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Abstract. Shock bound interaction zones (SBIZs) are ubiquitous in astro- 
physics. We present numerical results for 2D and 3D, plane-parallel, infinitely 
extended SBIZs. Isothermal settings and parameterized cooling are considered. 
We highlight and compare characteristic of such zones. We emphasize the mu- 
tual coupling between the turbulence within the SBIZ and the confining shocks, 
point out potential differences to 3D periodic box studies of supersonic turbu- 
lence, and contemplate on possible effects on the X-ray emission of such zones. 



1. Introduction 



Supersonically turbulent, shock-bound interaction zones (SBIZs) are impor- 
tant for a variety of astroph ysical objects. They contribute to structure for- 
mation in molecular clouds ( Hunter et al. _ 19861 : Audit and Hennebelld 120051 ; 



Vazquez- Semadeni et al.\ 120061; iHennebelle et al\ 120081 ') and to galaxy forma- 
tion ( Anninos and Nor man! 11996 ; Kang et al. 20051). They affect t he X-ray 
emission of hot-star winds (jQwocki et alj [19881 ; iFeldmeier et al. 1997) and the 



physics and emitted spectrum of colliding wind binaries ( Stevens et al. _ 19921; 



Nussbaumer and Walder 19931 ; Myasnikov and Zhekov 1998 : Folini and Walder 
2000). A large number of numerical and theoretical studies h ave great l y im- 
proved our understanding of SBIZs under various conditions (I Vishniad 1 1994; 

Blondin and Markdll996l : IWalder and Folinill 19981 : iHeitsch et aUl2005l : lFolini and Walder 
2006) 

In this paper, we want to highlight and compare some of the characteristics 
of such interaction zones under different conditions: isothermal - radiatively 
cooling, 2D - 3D, symmetric - asymmetric settings, early times - late times. 



2. Model Problem and Numerical Method 

We consider 2D and 3D, plane-parallel, infinitely extended SBIZs. Both, isother- 
mal and radiatively cooling settings are investigated. Two high Mach-number 
flows, oriented parallel (left flow, subscript I) and anti-parallel (right flow, sub- 
script r) to the x-direction, collide head on. The resulting SBIZ is oriented 
in the y-z-plane. It consists of a cold, dense layer (CDL) and, in the case of 
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Figure 1. Sketch of the SBIZ. Radiatively cooled material piles up in the 
CDL. Two hot post shock layers (LI and L2) exist in the case of parameterized 
cooling. Under isothermal conditions the SBIZ consists only of the CDL. 



radiative cooling, two layers of hot shocked gas (Fig. [[]). We investigated this 
system within the frame of Euler equations. In the polytropic equation of state 
we take 7 = 1.000001 (isothermal) or 7 = 5/3 (radiatively cooling). The right 
hand side of the energy equ ation we set to zero (iso thermal) or to a parameter- 
ized radiative loss function (TWalder and Folinilll996l ). Codes from the A-MAZE 
code package ( Walder and Folini 2000al ) are used. A coarse mesh is used for the 
upwind flows, a finer mesh for the CDL. Typical refinement factors are 2 8 in 
2D and 2 4 in 3D, resulting in 1280 (2D) and 256 (3D) cells in the y-direction. 
The meshes adapt automatically to the spatial extension of the CDL. The so- 
lution (density) on two adjacent refinement levels is compared, cells where the 
difference exceeds a prescribed tolerance are nagged. Refinement is applied to 
rectangular blocks that contain all flagged cells and typically also some unflagged 
cells. 



3. Characteristics of Interaction Zones 

Isothermal SBIZ: Average quantities, like slab thickness, auto-correlation 
length of the confining shocks, or energy loss per unit volume in the CDL, evolve 
approximately self-si milarly for symmetric flow collisions, where p\ = p r = p n 
and Mi = M r = M u (jFolini and Walderl 120061 ) . In 2D and 3D the root mean 



square Mach number M rms of the CDL scales linearly with M u , the mean density 
p m is independent of M u . This contrasts with the ID case, where p m oc M 2 . 

An intricate, 'self-regulating' interplay exists between the turbulence within 
the CDL and its driving: f e g = 1 — M{? ms with (3 ~ —0.6 in 2D. The fraction f e g 
of the upstream kinetic energy that passes the confining shocks unthermalized 
and M rms affect each other. The larger M u , the larger M rms , and the larger 
the fraction of the upstream kinetic energy that is thermalized only within the 
CD L. The dependence persis ts for asymmetric settings, where pi 7^ p r and M\ 7^ 
M r (|Folini and Walderl I2OO6I L In 3D we find the confining shocks to be steeper 



than in 2D (Fig. [21 left), thus f e g and M rms are larger while p m is smaller. 

The CDL is patchy. Patch sizes increase with the x-extension of the CDL 
and with decreasing M u , as does the auto-correlation length of the confining 
shock s. The density variance follows a 2 (ln{p)) = ln(l + M r 2 ms /4) (jPadoan et all 
19971 s ) up to M rms « 5 but then levels off (Fig. El right). This may be related 
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Figure 2. Left: Isothermal SBIZ in 3D for M u = 22. Shown is density, the 
confining shock is colored according to the local value of / e ff. On average, the 
confining shocks are steeper in 3D (/ G ff « 0.83 for M u = 22) than in 2D (f e g ss 
0.59 for M u = 22). Right: <j(ln(p)) versus M rms from 2D symmetric slab 
simulations ( red diamonds) . The 2D slab resul ts deviate from published fitting 
functi ons bv lPadoan et al\ (|1997| ) (black) and lPassot and Vazquez- Semadenil 
(|1998f) (green) for large M rms . 




to the strongly anisotropic velocity field within the CDL, with larger (smaller) 
velocities along (perpendicular to) the upstream flow direction (Fig. El left). 



One sided cooling SBIZ: Astrophysical examples include wind blown bub- 
bles or planetary nebula, where a fast wind collides with a slow precursor, form- 
ing a hot (> 10 7 K) reverse and a cool (< 10 6 K) forward shock. A charac- 
teristic feature of SBIZs where one of the confining shocks cools efficiently but 
not the other are clumps of cold material that b reak off from the CDL a nd drift 
into the non-cooling, hot post-shock material (jWalder and Folrnil Il998l ). Den- 
sity images show a striking similarity with HST pictures of the Helix nebula. 
A prerequisite is that the corresponding boundary of the CDL is occasionally 
hit by a shock wave, which triggers a Richtmeyer-Meshkov instability. The 
shock wave can be generated by a thermal instability of the oth er hot post 
shock layer ( Strickland and Blondin 1995I ; IWalder and Folini 19961 ) or by inho- 
mogeneities in the upstream flow. The CDL is weakly turbulent, the velocity 
field is isotropic. 



Two sided cooling SBIZ: If both confining shocks cool efficiently with spa- 
tially resolved hot post shock layers (cooling limit = upstream temperature), 
the evolution of the CDL at early times can become much more violent than in 
the isothermal case (Fig. HI left). At later times, the hot post shock zones rather 
exert a cushioning effect (Fig. HI middle and right). The characteristics of the 
CDL are hardly affected. The velocity pdfs show not much change from early 
to late times (Fig. O right). Much more obvious is the difference to isothermal 
simulations, where the velocity field within the CDL is much more anisotropic 
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Figure 3. Velocity pdfs (x- and y-direction, two times) for an isothermal 
(left) and a radiatively cooling (right) 2D simulation with M u = 22 or v u — 
1.8 ■ 10 7 cm/s. An early (bluish colors) and a late time (redish colors) are 
shown. Wider (narrower) pdfs represent x- velocities (y- velocities). Density 
plots of the cooling case are shown in Fig. [4] (left and right frame). 



(Fig. [3l left). Also, for the same upstream parameters, the root mean square 
Mach nu mber in the CDL is much lower in the cooling case than in the isother- 
mal case (Walder and Foliriill2000b : Folini and Warden l2006l ). Resolved cooling 
layers appear to have a damping effect on the turbulence within the CDL. 



4. Discussion and Conclusions 

Shock bound interaction zones come in many physical varieties and an intri- 
cate interplay exists between the CDL and the confining shocks. Much work 
remains to be done, especially also with regard to observational and theoretical 
consequences. The presented results raise, in particular, the following questions. 

One question concerns the wide spread assumption that the relative speed of 
high velocity colliding flows translates more or less directly into the hardness of 
the X-ray spectrum. At least under isothermal conditions the confining shocks 
are not perpendicular to the colliding flows but steepen with increasing M u . 
This should soften the spectrum and shift part of the em ission from the confining 
shocks to internal shocks of the CDL (jSmith et aZ.ll2OQ0l ) . Whether geometrically 
thin, resolved cooling layers behave similarly remains to be studied. 

Another topic coming to mind are the relation between SBIZs and 3D peri- 
odic box simulations. The later have enormously improved our understanding of 



supersonic turbulence (IMac Lowlll 999; Bo ldvrey et al. 2002; Pado an et a/.ll2004l ; 
Kritsuk and Normanl |2004| ; ISchmidt et all 120091 ; iFederrath et aZ.ll2009l ). All the 



more one may ask to what degree these results carry over to the CDL in a SBIZ, 
where surface effects play a role. Three points may be raised. 

In an isothermal SBIZ, M rms , the energy input into the CDL, and the spatial 
scale on which this energy input is modulated all depend on each other. In 3D 
box simulations, the energy input and its spatial modulation - the driving wave 
length - are usually chosen independently. The effect of this difference on the 
turbulence characteristics remains to be clarified. A second point is the velocity 
field, which is clearly anisotropic in the CDL but usually isotropic by design in 3D 
box simulations. Finally, when computing a CDL-averaged structure function 
(SF) for our 3D CDL (compute SFs for each point in the CDL, using only points 
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in the CDL, then taking the average over all SFs ) it resembles more closely the 
observatio n based SFs by 



3D boxes (jBoldvrev et al 



Gustafsson et al 



t han those typically found fo r 
120021 : iPadoan et al.ll200J : iKritsuk and Normadl2004h . 



The reason is not yet clear. It could be the finite spatial extension of the CDL, 
its inhomogeneous, anisotropic interior, or simply a too coarse resolution. 
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Figure 4. Two sided cooling SBIZ, time evolution, logarithm of density. 
CDL in red and white, hot post shock layers in light blue. Filamentary and 
dynamic at early times, the CDL becomes more 'bulky' later on. Nevertheless, 
small scale features at the surface of the CDL continue to form and vanish. 
No CDL was present at time 0. 

Further studies of SBIZs are needed, both with basic physics and with 
more elaborate models, and with particular attention to the interplay between 
the confining shocks and the turbulent CDL. 
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